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We present here the results of our relativistic many-body calculations of various properties of the 
first six low-lying excited states of indium. The calculations were performed using the relativistic 
coupled-cluster method in the framework of the singles, doubles and partial triples approximation. 
We obtain a large lifetime ~ 10s for the [4p®]5s^5p3/2 state, which had not been known earlier. Our 
results could be used to shed light on the reliability of the lifetime measurements of the excited 
states of atomic indium that we have considered in the present work. 
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I. INTRODUCTION 

Indium (In) has been laser cooled and trapped a few 
years ago Following this experiment, a proposal has 
been made to search for the permanent electric dipole 
moment (EDM) in this atom 0] . It would ineed be desir- 
able to carry out high precision measurements and many- 
body calculations of other properties of this atom. A few 
measurements of the magnetic dipole hyperfine structure 
constants of the first three low-lying states of In are al- 
ready available [1, However, the reported theoretical 
results obtained using different variants of the relativis- 
tic coupled-cluster (RCC) method at the singles, dou- 
bles and important triples excitations level (CCSD(T) 
method) are not able to reproduce them to within one 
percent accuracy 0, Q . This suggests that the role of 
correlation effects for this property is of crucial impor- 
tance. In addition, it would also be worthwhile to cal- 
culate different transition amplitues in In for a number 
of reasons. First of all, the behavior of the correlation 
effects in these properties could be quite different than 
in the hyperfine structure constants. Furthermore, these 
amplitudes in conjunction with the hyperfine constants 
can be employed to verify the accuracy of the wave func- 
tions for the proposed EDM calculations 0] or perhaps 
for the parity nonconservation in this system if at all 
an experiment to observe this effect can be carried out 
in this atom and also to determine the polarizabilities, 
lifetimes, oscillator strengths, branching ratios etc. for 
various states. 



In this work, we calculate the excitation energies (EEs) 
and different transition amplitudes due to allowed and 
forbidden electromagnetic transitions among the first six 
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FIG. 1: Schematic low-lying energy level diagrams and decay 
channels of the low- lying states in In. 



low-lying states, giving a total of 34 possible transitions 
(see Fig. [1]), using the relativistic CCSD(T) method. 
These results are further used to determine transition 
rates, branching ratios and lifetimes of the above states. 
These properties are also important from an astrophys- 
ical point of view 0| ■ Safronova et at. have reported 
EEs and electric dipole (El) transition amplitudes for 
a number of states and compared with previous calcu- 
lations and measurements [5,]. However, they have only 
considered these transition amplitudes to estimate the 
lifetimes of various states, but contributions from the for- 
bidden transitions are not included. In our calculations, 
we have taken into account the forbidden transition am- 
plitudues in the evaluation of the lifetimes of different 
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states. 



II. THEORY AND METHOD OF 
CALCULATIONS 

The transition rates (in s~^) due to various transitions 
are given by @ 



and 



AMI 



aB2 



iM2 



2.02613 X IQis 
A3(2J/ + 1) 

2.69735 X 10" 
A3 (2 J/ + 1) 

1.11995 X 10^8 
A^(2J/ + 1) 

1.491 X 10" 



cEl 



rtMl 



cE2 



5 



M2 



A5(2J/ + 1) f^'' 



(2.1) 
(2.2) 
(2.3) 

(2.4) 



where A (in A) and Sf^^{^ l(/l|0||i)P) (in atomic unit 
(au)) are the wavelengths and hne strengths due to the 



corresponding transition operator O, respectively. 

The lifetime (r/) of a given state / is just the reciprocal 
of the total transition rate of that state due to all possible 
transition channels; i.e. 
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where A'^^ - is the transition rate due to operator O and 
sum over i and O represents the total transition rate from 
state / to all possible states i and due to all possible 
operators. 

The branching ratios due to an operator O from a state 
/ due to the lower states are given by 



(2.6) 



To evaluate the line strengths, we use the following 
reduced matrix elements at the single particle orbitals 
level for the El, Ml, E2 and M2 operators [i| 





l|el| 






ml 1 


1 




||e2| 


1 i^t) 


and 






Wll 


m2 1 


1 i^t) 



{nj II II n,) drriiPfP, + QfQ,) - 

{~^f \\C<'^ \\k,) r drr^-^l±^{PfQ, + QfP,), 
Jo ^ 



(«/ II C^^'^ II '^^) / drr-'iiPfP, + QfQ,) - — 
Jo 'Q: 



-{PfQr + QfPi) + {PfQ. - QfP^)]), (2.7) 

(2.8) 

— {PfQ. + QfP^) + {PfQ^ - Q/P.)]}(2.9) 



3a 



iPfQ^ + QfP^), 
I 



(2.10) 



where, j's and k's are the orbital and relativistic angu- 
lar momentum quantum numbers, respectively, P and Q 
represent the radial parts of large and small components 
of single particle Dirac orbitals, respectively, lj — €f ~ €i 
for the orbital energies es, a is the fine structure constant 
and the reduced Racah coefficients are given by 

wiicwii«,) = (-i)^/+vy(2j/ + i)(2j, + i) 



with 



7r(Z, TO, I') 



1 for ? + TO + V 
otherwise. 



(2.12) 



In the above expressions and in the remaining part of 
the paper, we have used au unless they are mentioned 
explicitly. 



In order to determine the above properties, we cal- 
culate the atomic wave function (j^'t,)) with a valence 
orbital v by expressing it in the RCC framework as 



(2.13) 



where we define a reference state 1$^,) by appending the 
appropriate valence orbital v to the Dirac-Fock (DF) 
wave function d'l'o}) with the configuration similar to 
cadmium; i.e. [Ap^]Ad^^bs'^ . Here T and represent the 
excitation operators due to core-core and core-valence 
electron correlations. In the CCSD(T) method, the T 
and operators are defined as 



T = Ti+T2 and 5*^, = Si^ + S2V, 



(2.14) 



where 1 and 2 in the subscripts represent for single and 
double excitations, respectively. 

The equations determining the coupled-cluster ampli- 
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tudes and energy can be expressed in compact forms as 

($^|{^}|$0) = 6o,LAE,orr (2.15) 

and 

{<f^\{He^}{l + S^}\^,,) = ($^|1 + 5,|$„) 

($„|{i??^}{l + S„}|<i>,) 
= ($^'|<5l^.„ + 5,|$,)AS,,(2.16) 

where the superscript L{= 1,2) represents for the excited 

hole-paerticle states, He^ denotes the connected terms 
of the Dirac-Coulomb (DC) Hamiltonian with the T op- 
erators, AEcorr and AEy are the correlation energy and 



attachment energy (also equivalent to negative of the ion- 
ization potential (IP)) of the electron of orbital v, respec- 
tively. The reference states in Eq. ()2.15p and Eq. (j2.16p 
contain different number of particles, hence the Hamil- 
tonian used in the respective equations describe differ- 
ent number of particles in our Fock space representa- 
tion. Contributions from the important valence triple ex- 
citations are included perturbatively through the above 
equations. 



The transition matrix element of a physical operator 
O between l^*/) and j^'i) in our approach is given by 



(*/|0|*.) 



{'i>f\{l + S}}0{l + S,}m 
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where O = e^^ Oe^ and N = e^'' are two non- 
truncating series in the above expression. We evaluate 
them by considering terms whose leading contributions 
arise in fourth order perturbation theory or lower. Con- 
tributions from to the normalization of the wave func- 
tions {N) are determined explicitly as follows 
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We have used Gaussian type orbitals (GTOs) to con- 
struct the single particle orbitals for the Dirac-Fock 
(|$o)) wave function. The large and small components 
of the Dirac orbitals in this case are expressed as 



and 



(2.19) 



(2.20) 



where the summation over k is for total number of GTOs 
used in each symmetry, and are the normalization 
constants for the large and small components, respec- 
tively, and we use (^ + 7) operator to expand the small 
component Dirac orbitals to maintain the kinetic balance 
condition with its large component. In the present cal- 
culations, we have considered 9 relativistic symmetries 
(up to g symmetry) and 28 GTOs for each symmetry to 
generate the orbitals. In order to optimise the exponents 



to describe orbitals from various symmetries in a smooth 
manner, we use the even tempering condition 



ak = aol3' 



fe-i 



(2.21) 



where ao and /? are two arbitary parameters that can 
be chosen suitably for different symmetries. We have 
considered ao = 7.5 x 10"'* for all the symmetries and 
/3 are taken as 2.53, 2.45, 2.58, 2.75 and 2.83 for s, p, d, 
f and g orbitals, respectively. For the RCC calculations, 
we have considered excitations up to first 16s, 16p, 16d, 
14/ and 13g orbitals as the remaining orbitals have large 
continuum energies. 



III. RESULTS AND DISCUSSION 

In accordance with Koopman's theorem, the energies 
of the virtual orbitals obtained in our calculations are the 
IPs at the DF level, since our DF wave function is com- 
puted using the closed-shell configuration [4p^]4d^°5s^. 
We present the IPs from NIST 10], from other calcula- 
tion as well as our calculations in Table HI 

As observed from Table HI our IP results are within 
0.5% except for the 5d states (which are around 1% ac- 
curate) compared with the results given in (lOj . In an 
earlier work, Safronova et al. have reported the results 
for these quantities based on the linearlized version of 
the relativistic CCSD (T) method using a B-spline ba- 
sis 0. The major differences between this and our work 
are the different basis sets used in the two calculations 
and the additional non linear clusters in our calculation. 
Our results are in better agreement with the high preci- 
sion NIST results than those of 0] for all the states that 
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TABLE I: IPs (AEys) of different states of In in cm"\ Ab- 
solute error of our CCSD(T) results compared to the quoted 
results in are given as A. 



State 


NIST 
in cm~^ 


''Others 
in cm~^ 


"^Koopman 
in cm~^ 


CCSD(T) 
in cm~^ 


A 

in % 


5pi/2 


46670.11 


46189 


41521.74 


46581.47 


0.19 


5P3/2 


44457.51 


44031 


39522.20 


44361.04 


0.22 


6Sl/2 


22297.15 


22442 


20567.70 


22291.74 


0.02 


6pi/2 


14853.21 


14833 


13977.92 


14819.07 


0.23 


6p3/2 


14554.89 


14532 


13718.04 


14519.46 


0.24 


5^3/2 


13777.90 


13581 


12389.68 


13633.48 


1.05 


5^5/2 


13754.57 


13554 


12373.64 


13603.88 


1.10 



" Reference [iG]. 
^ Reference [g]. 
" This work. 



we have considered. Given the high accuracy of our IPs 
and therefore the excitation energies, we can accurately 
determine the wavelengths for various transitions in or- 
der to determine the transition rates and the lifetimes of 
different excited states. We can also use the wavelengths 
from NIST data to obtain the lifetimes and compare them 
with the results from the relativistic CCSD(T) method. 

In Table [TTl we present the line strengths obtained us- 
ing the DF and relativistic CCSD(T) methods for both 
the allowed and fobidden transitions. Safronova et al. 
have given the results only for the allowed transitions 
and they have not verified explicitly the contribu- 
tions from the forbidden transitions. In fact, our tran- 
sition strengths for the allowed transitions differ slightly 
from theirs and the cause of the differences between the 
approximations employed in the two cases have already 
been discussed earlier. The line strengths due to the for- 
bidden transitions are not small in many of the cases. 
The influence of the forbidden transitions should be ver- 
ified in the determination of transition rates, BRs and 
lifetime estimations as they may be important in some 
cases. 

We present the wavelengths, transition rates, branch- 
ing ratios and lifetimes of different states of In in Table 
mil These quantities are determined using both the cal- 
culated wavelengths that are estimated from the excita- 
tion energies obtained in this work and the experimental 
wavelengths from NIST data [l0|. The ab initio results 
are given as I and wherever the experimental wavelengths 
used are given as II. We also give measured lifetimes re- 
sults based on different experimental techniques jT2l - [l7j 
in the same table. 

The difference between the experimental results and 
obtained calculations for EEs are treated as possible un- 
certainties associated with them, which are given in per- 
centage in Table m Uncertainties in the calculated transi- 
tion matrix elements are obtained by finding out the con- 
tributions from higher angular momentum orbitals using 



TABLE II: Line strengths (in an) due to allowed and forbid- 
den transitions between different states in In. Numbers given 
in the parentheses and square brackets represent estimated 
errors and powers in 10, respectively. 



Transition DF CCSD(T) Others [5] 



the dominant many-body perturbation diagrams and are 
mentioned in the parenthesis of the results presented in 
Table im In the final lifetime estimation of various states, 
we consider central values given as II in Table IIIII and 
uncertainties are determined from the above error bars. 
These results are reported as recommended values (Reco) 
in Table Uni 

As we find from the above table, the branching ratios 



5^5/2 


El^ 


6P3/2 


251.95 


188(2) 


186 


5c?5/2 


M2^ 


■ QV'i /o 

^fo 1 2 


4.9f3l 


3.8fl)[3l 




5c?5/2 


M2^ 


■ 6pi/2 


893.6 


685(5) 




5c?5/2 


El^ 


5z>^ /o 


20.55 


16.3(5) 


15.2 


5^5/2 


M2^ 


■ 5p3y2 


399.59 


474(10) 




5^^5/2 


M2^ 


' 5jDi /o 


62.44 


92(3) 




5c/r Jo 
-'^o / 2 


Ml^ 


'-'^6 1 2 


2.40 






5c?5/2 


E2^ 


56?3/2 


5.0[3] 


2.9(2)f3l 




5(^5/2 


E2^ 


6Sl/2 


9.4[3l 

L J 


6.6(l)[3l 

V / L I 




5^3/2 


El^ 


6P3/2 


27.88 


20(1) 


20.5 


5^3/2 


M2 


■ 6p3/2 


0.0 


3.5(l)[-3l 

V / L J 




5c^3/2 


El^ 


6pi/2 


139.07 


104(5) 


103 


5^3/2 


M2^ 


• 6pi/2 


83.44 


65(3) 




5c?3/2 


El^ 


5?>3/2 


2.30 


1.84(2) 


1.71 


5c?3/2 


M2^ 


■ ^P3/2 


0.0 


0.27(1) 




5(^3/2 


El^ 


5pi/2 


9.84 


7.7(4) 


7.24 


5(^3/2 


M2^ 


' 5zPi /o 


5.90 


4.9(3) 




5*^3/2 


Ml^ 


■ 6Si/2 


4.8[-12l 

L J 


2(l)h9l 




5c?3/2 


E2^ 


6-51/2 


6.3[3] 


4.4(l)[3l 




6P3/2 


Ml^ 


■ 6pi/2 


1.33 


1.33(1) 




6p3/2 


E2^ 


6pi/2 


1.64f4l 

L J 


1.32(1) f4l 

V / L J 




6p3/2 


El^ 


6Sl/2 


88.96 


72.9(1) 


70.3 


6p3/2 


Ml^ 


■ 5p3/2 


4.9f-9l 


2.1(5)[-4l 

\^/ L J 




6p3/2 


E2^ 


5p3/2 


131.93 


106(8) 




6p3/2 


Ml^ 


■ 5pi/2 


9.2[-4] 


6(l)[-4] 




6p3/2 


E2^ 


5pi/2 


96.90 


77.3(8) 




6pi/2 


El^ 


6Sl/2 


45.81 


37.5(1) 


36.1 


6pi/2 


Ml 


• 5p3/2 


1.0[-3] 


1.4(3)[-3] 




6pi/2 


E2^ 


5P3/2 


149.33 


120(7) 




6pi/2 


Ml 


• 5pi/2 


3.9[-10] 


1.2(8)h5l 




6Si/2 


El^ 


5P3/2 


11.26 


8.8(2) 


8.56 


6Sl/2 


El^ 


5pi/2 


4.68 


3.67(2) 


3.64 


5p3/2 


Ml^ 


• 5pi/2 


1.33 


1.31(1) 




5p3/2 


E2^ 


5pi/2 


236.42 


181(1) 
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TABLE III: Wavelengths (A) in A, transition rates (A) in , branching ratios (F) and lifetimes (r) in ns for the considered 
excited states in In. We consider the calculated and experimental values of As to determine the above quantities which are 
given as I and II, respectively. We present the recommended (Reco) values for rs after accounting possible errors and compared 
them with their experimental (Expt) results. Numbers given in the parentheses and square brackets represent estimated errors 
and powers in 10, respectively. 



Upper 
state (/) 


Lower 
state (i) 


Channel 



A 

I 


II 


I II 




I 


II 


Reco 


Expt 


5^5/2 


6P3/2 


El 


1.09[5] 


1.25[5] 


48781.17 


32544.0 


0.0002 


6.22 


6.27 


6.3(2) 


"7.6(5) 




6p3/2 


M2 








1 [—1 nl 


- 0.0 








''7.1(6) 




6pi/2 


M2 


8. 23 [4] 


9.10[4] 


4 5f— 51 


2 7f— 51 


- 0.0 












5P3/2 


El 


3251 


3257 


1.61 [8] 


1.60 [8] 


0.9998 












5p3/2 


M2 






3.2[— 3] 


3.2[— 3] 


- 0.0 












5pi/2 


M2 


3032 


3038 


8.9-4 


8.8-4 


- 0.0 












5^3/2 


Ml 


3.38[6] 


4. 29 [6] 


2.8[-7l 


1.4[-7l 


~ 0.0 












5^3/2 


E2 






1.2[-12] 


3.7[-12] 


- 0.0 












6Si/2 


E2 


1.15 [4] 


1.14[4] 


6.10 


5.61 


- 0.0 










5d3/2 


6P3/2 


El 


1.13[5] 


1.29 [5] 


7316.78 


4935.01 


3.0[-5] 


5.99 


6.03 


6.0(3) 


"6.3(5) 




6p3/2 


M2 






7 If— 161 




- 0.0 








^^7.0(4) 




6pi/2 


El 


8.43[4] 


9. 30 [4] 


8. 78 [4] 


6.55[4] 


4.0[-4] 












6pi/2 


M2 




5.7[-ll] 


3.5[-ll] 


- 0.0 














5P3/2 


El 


3254 


3259 


2.71 [71 

' L ' J 


2. 69 [7] 


0.16 












5P3/2 


M2 






2.8f-6l 


2.7f-6l 


- 0.0 












5pi/2 


El 


3035 


3040 


1.40 [8] 


1.39 [8] 


0.84 












5pi/2 


M2 




7.1[-5] 


7.1h5l 

L J 


~ 0.0 














6Si/2 


Ml 


1.16 [4] 


1.17[4] 


l.l[-8] 


1.0[-8] 


~ 0.0 












6Si/2 


E2 






5.99 


5.53 


~ 0.0 










6P3/2 


6pi/2 


Ml 


3.34[5] 


3. 35 [5] 


2.4[-4l 

L J 


2.4[-4l 

L J 


~ 0.0 


57.67 


58.34 


58(1) 


''55(4) 




6pi/2 


E2 






8.9[-7] 


8.7[-7] 


~ 0.0 












6Si/2 


El 


1.29 [4] 


1.29 [4] 


1.73 [7] 


1.71[7] 


0.9999 












5P3/2 


Ml 


3351 


3344 


0.04 


0.04 


~ 0.0 












5p3/2 


E2 






69.99 


70.73 


4.1[-6] 












5pi/2 


Ml 


3119 


3114 


0.13 


0.13 


- 0.0 












5pi/2 


E2 






73.34 


73.93 


4.3[-6] 










6pi/2 


6Si/2 


El 


1.34 [4] 


1.34 [4] 


1.58[7] 


1.56[7] 


0.9999 


63.16 


63.81 


63.8(8) 


'*55(4) 




5P3/2 


Ml 


3385 


3378 


0.49 


0.49 


- 0.0 












5P3/2 


E2 






150.79 


152.35 


9.7[-6] 












5pi/2 


Ml 


3148 


3143 


0.005 


0.005 


~ 0.0 










6Si/2 


5p3/2 


El 


4531 


4513 


9.58[7] 


9. 70 [7] 


0.6431 


6.71 


6.63 


6.6(2) 


"7.5(7) 




5pi/2 


El 


4117 


4103 


5. 33 [7] 


5. 38 [7] 


0.3569 








=7.0(3) 
























^ 7.4(3) 


5P3/2 


5pi/2 


Ml 


4.50[4] 


4.52[4] 


0.0967 


0.0957 


0.9969 


10.31[9] 


10.42[9] 


10.4(2) [9] 




5pi/2 


E2 






0.0003 


0.0003 


0.0031 











"Reference 
''Reference 
'^Reference 
''Reference 
Reference 
^Reference 



due to the allowed transition channels completely dom- 
inate over the forbidden transition channels. Therefore, 
the lifetimes of the excited states except for the 6^3/2 
state are almost entirely determined by the allowed tran- 
sitions. The 5p3/2 state is the fine structure partner of 
the the ground state, its lifetime is determined from the 
forbidden transitions. For this case, the Ml transition 
clearly dominates over the E2 transition and that is evi- 



dent from their branching ratios. We obtain a large life- 
time, ^ 10s, for this state. The lifetime of the 5^5/2 state 
obtained from our calculation is in reasonable agreement 
with the available experimental data. We find that the 
lifetime of this state is almost entirely due to the El de- 
cay channel to 5p3/2 state. As the wavelength of the 
transition 5c?5/2 6^3/2 is very large, the branching ra- 
tio of the 5(1^/2 state is small. Our calculated lifetime 
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for the 5^3/2 state agrees well with the experimental re- 
sults. We find 84% and 16% branching ratios from this 
state to the ground and 5j33/2 states, respectively through 
the El channel. Contributions from the forbidden tran- 
sitions are also negligible in this case. Similar agreement 
between our calculated and experimental results for the 
lifetimes of the 6p states are found, but the experimental 
results have large error bars compared to our calcula- 
tions. There is also a marginal difference between the 
measured lifetimes and the calculated lifetimes of the 6s 
state, although they are within the common error bar. 
The branching ratios from this state to the ground and 
5p3/2 states are of the order of 35% and 64%, respectively. 
This trend is different for the 5c?3/2 state as discussed 
above. 



IV. CONCLUSION 

We have estimated the branching ratios and life- 
times of certain low-lying excited states of indium. We 



have carried out calculations of the excitation energies 
and line strengths using the relativistic coupled cluster 
method. We have also compared our ab initio results 
with the results obtained using the experimental wave- 
lengths and measured lifetimes from different experimen- 
tal techniques. We find that the forbidden transitions do 
not contribute significantly to the lifetimes of most of the 
states that we have considered. A large lifetime for the 
5P3/2 state 10s) has been found from this work, which 
is completely due to the forbidden transitions to its fine 
structure partner; the ground state. 
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